Domain walls in gapped graphene 
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The electronic properties of a particular class of domain walls in gapped graphene are investigated. 
We show that they can support mid-gap states which are localized in the vicinity of the domain 
wall and propagate along its length. With a finite density of domain walls, these states can alter 
the electronic properties of gapped graphene significantly. If the mid-gap band is partially filled,the 
domain wall can behave like a one-dimensional metal embedded in a semi-conductor, and could 
potentially be used as a single-channel quantum wire. 



Graphene is a one-atom thick layer of carbon atoms 
with a hexagonal lattice structure and where electrons 
within ~ lev of the Fermi energy obey a Dirac equa- 
tion and have a linear dispersion relation u = with 
Fermi velocity vf ~ c/300. It has been studied as an 
analog of relativistic field theory [l[ where relativistic 
quantum mechanics and field theory phenomena special 
to 2+1 space-time dimensions 00 could be realized in 
nature. It was identified and studied in the laboratory in 
2004 0] and it turns out to have high conductivity and 
carrier mobility and other interesting properties which 
make it apromising material for applications in electronic 
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devices 

An important current problem is to modify graphene 
so that it has a gap in its energy spectrum Q. The idea 
is that, with a small gap, so that electrons obey a Dirac 
equation with mass, it would retain the good features of 
graphene and could also be used as a semi-conductor in 
applications where a gap is essential, for example, a field 
effect transistor. Breaking the symmetry which inter- 
changes the two triangular sublattices of the hexagonal 
graphene lattice will gap the spectrum [l[ . This could be 
accomplished, for example, by giving electrons residing 
on A sites a different energy from those on B sites by in- 
troducing a staggered chemical potential [l[. It can also 
arise from deformations of bonds on the graphene lattice 
Q analogous to those known from the study of carbon 
nanotubes [ljj. A third possibility is to use multi-layer 
graphene where the layers can be stacked so that their in- 
teraction breaks the sublattice symmetry. In all cases, to 
retain the features of the Dirac equation the gap should 
be much less than the nearest neighbor hopping ampli- 
tude t ~ 2.7ev. 

The diatomic material Boron-Nitride (BN) has the 
same lattice structure and valence electrons as graphene 
and a staggered chemical potential by virtue of having 
different atoms on the two sublattices. Monolayers have 
been made in the laboratory [ll[. However, the gap is 
too large ~ A.bev for Dirac electrons. An approach cur- 
rently being pursued is to attach a graphene monolayer 
to a BN substrate. The resulting gap in graphene is es- 
timated to be ~ 53mev [l2| which is in the interesting 
range. Another approach is epitaxial growth of graphene 
on a Silicon-Carbide substrate where a larger magnitude 
gap ~ 2.6ev has been observed [HI]. 



In this Paper, we shall consider line-like domain wall 
defects in the mass pattern in graphene which is gapped 
by a sublattice symmetry breaking staggered chemical 
potential. We shall find that they can have significant 
electronic properties. The domain walls are shown in 
Figs. 1(b), 1(c) where the zig-zag and armchair walls form 
boundaries between regions where the staggered chemical 
potential is shifted between the two sublattices. Such do- 
main walls could be realized naturally in BN and would 
be inherited by graphene on a BN substrate, for example. 
We shall show that they can give rise to a band of mid- 
gap states. These states are localized in the vicinity of 
the wall and propagate along its length. If the mid-gap 
band is partially filled, the domain wall can behave like 
a one-dimensional metal embedded in a semi-conductor, 
and could potentially be used as a single-channel quan- 
tum wire. One might imagine that, once techniques for 
deposition of graphene monolayers on substrates are bet- 
ter developed, the conditions for existence of these do- 
main wall wires could be created and manipulated to the 
point where they could be used to print electric circuits 
on graphene sheets. 

Mid-gap states already play an important role in 
graphene. It was pointed out long ago [l[ that an in- 
dex theorem governs the degeneracy of the E = Lan- 
dau level in graphene in a magnetic field and this level is 
half-filled in the neutral material. This observation has 
spectacular experimental confirmation in the half-odd- 
integer anomalous quantum Hall effect [f| In ad- 
dition, theoretical studies of point-like vortex defects in a 
mass condensate due to a Kekule distortion of graphene 
find mid-gap electron states which can give the vortices 
fractional charge [lj][l6|-[l!|, thus giving a two dimen- 
sional realization of a phenomenon previously known to 
occur in one-dimensional linearly conjugated polymers 
such as polyacetylene [2(|-[22]. Similar states bound to 
vortices in a proximity-induced superconducting conden- 
sate in graphene could lead to anyonic statistics with 
potential applications to quantum computing [23| . An 
essential common feature of these examples is the exis- 
tence of "zero-mode" mid-gap states in the spectrum of 
the Dirac Hamiltonian which arise from the interaction 
with fields that have a non-trivial topology. In the case 
of the vortex, this topology is due to the vorticity. Let 
us consider a simple example to show that a related phc- 
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nomenon takes place for a domain- wall. Consider the 
4 x 4-matrix graphcnc Dirac Hamiltonian with the addi- 
tion of a mass term (which might arise from a staggered 
chemical potential): 



H = hvp 
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The two diagonal blocks correspond to the two graphcnc 
valleys, which transform into each other under parity and 
time reversal. A domain wall is described by replacing 
the mass m in Eq. |T]) by a function m(x) which depends 
on one of the coordinates, x with a soliton profile 



lim m(x) 



-m < 



lim m(x) = m > (2) 



The energy spectrum then has the same gapped conduc- 
tion and valence band branches as would occur if m(x) 
in Eq. (fTJ) were a constant with the asymptotic value of 



the mass, 
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These describe 



electrons in the bulk semi-conductor away from the wall. 
As well, there is a gapless mid-gap branch whose wave- 
functions have support near the wall. Explicitly, the (un- 
normalized) wave-functions and eigenvalues are 



/ / \ ik^y/h—^r ff dx'm(x') 

lpR(X, y) — e vUI ~ J o *• ' 



E = VFky (3) 



E = -v F k y (4) 



Note that this solution exists and is continuum normal- 
izable for whatever the profile of the position-dependent 
mass term, it only needs to have the asy mptotic behavior 
of a topological soliton as in Eq. ^ [24(. In particular, it 
should be applicable to a one lattice spacing thick domain 
wall such as those drawn in Figs. 1(b), 1(c). 

What we have found are two bands of mid-gap states 
corresponding to one left- and one right-moving one- 
dimensional massless fermion (for each spin degree of 
freedom) traveling along the length of the domain wall. 
An effective Lagrangian describing them would be 



s (d t - v F d x )tp Ls + 4>R S {dt + v F d x )ipR S (5) 



where s labels the two spin states. Effects of impuri- 
ties and local interactions can be important in one di- 
mension and should be taken into account. Four- fermion 
operators are perturbatively marginal and adding those 
which do not implement umklapp processes yields the 
Tomonaga-Luttingcr model which is a solvable confor- 
mal field theory with well-known properties. 
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FIG. 1: a) A hexagonal graphene lattice with triangular sub- 
lattices A (black dots) and B (white dots) connected by vec- 
tors si = (0, -l)a,S2 = (v / 3/2, 1/2) a, s 3 = (-V^/2,1/2) a 
with lattice constant a = 1.23A. b) A zig-zag domain wall, 
c) An armchair domain wall. In b) and c), the chemical po- 
tentials are — /i on black dots and +/i on white dots; the 
sublattices which have higher/lower chemical potentials are 
interchanged at a domain wall creating a line of miss-matched 
neighbors, denoted by a dashed line. 



To understand the structure of the bands in more de- 
tail, we must take a closer look at the tight-binding lattice 
model of gapped graphene. We shall find features which 
are not reflected in the continuum analysis which is only 
valid in a small region near E = 0. The Hamiltonian is 



B,i 

+ E ^ a A a A - E ^ b B b B 



A,i 



(6) 



where a A ,aA,b^ B , bs are the quantum amplitudes for an 
electron to occupy sites labeled A and B on the sublat- 
tices A and B, respectively. The lattice and sublattices 
are depicted in Fig. 1(a). The first terms in Eq. ([6]) de- 
scribe electron tunneling between nearest neighbor sites. 
The terms proportional to fj, are on-site energies. They 
break the sublattice symmetry and generate a gap. The 
Schrodingcr equation is 

(E - n)a A = t J2 b A+s, , (E + /-i)b B = t OB-*, (7) 
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To study the zig-zag domain wall in Fig. 1(b), we solve 
Eq. ([7]) with fj, replaced by /i sign(A y ) and (i sign (B y ). 
The spectrum has branches corresponding to electrons 
propagating in the bulk of the gapped graphene away 
from the wall, 
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2 

(8) 

Here, {k x , k y ) are wave-vectors. This bulk spectrum has 
a gap 2/i and is symmetric about E = 0. Then, there are 
two branches with wave-functions which fall off exponen- 
tially with transverse distance |^4 y | , | i?^ | from the wall 
and are oscillating functions of the longitudinal A X ,B X - 
coordinates with wave- vector k v : 



E 



E = t- 




-2 V ^ a t 



cos^—hx < (9) 



4*2 C0S 2 y^ kx ; cos ^: kx > 0(10) 



The band in Eq. (J9]) is located inside (and reaches slightly 
below) the negative energy bulk states Eq. ([5]) . The other 
branch Eq. (fTC)|) covers the interval [t — \J [i 2 + At 2 , t — fj] . 
If t > /i, this band crosses zero energy (E = 0) at two 
values of k x , and thus agrees with the continuum analysis 
which is only valid when /i << t and which predicts the 
existence of two zero energy modes - one for each cross- 
ing. The spectrum and density of states for fi = O.bt 
are depicted in Fig. 2. Note that, unlike the spectrum of 
bulk graphene Eq. ([5]) the zig-zag domain wall spectrum 
is not symmetric about E = 0. This is evident from its 
structure displayed in Fig. 1(b), where the miss-matched 
sites along the wall are entirely black dots with chemi- 
cal potential — fi. The zig-zag domain wall violates the 
symmetry which reflects the sign of the energy There is 
an anti-wall where the miss-matched bonds are entirely 
white dots - with energy +/i. Its domain wall spectrum 
would have opposite sign to Eqs. © and (JTUJ) ■ 

We can get an intuitive understanding of the spectrum 
in Eqs. ([9]) and (fT0|) in the limit where /j, is large. Initially, 
neglecting t, there are two energy levels, /i for an elec- 
tron sitting on a white dot and —\i for an electron sitting 
on a black dotin Fig. 1(b). Then, if we turn on small t, 
the largest effect is for the black dots on each side of the 
domain wall which have a nearest neighbor at the same 
zcroth order energy, —fx. Turning on the hopping would 
split the degeneracy of these sites to — \x + t and — fj, — t. 
Note that this does not happen for sites in the bulk away 
from the domain wall, since they are not degenerate with 
their neighbors - corrections to their spectrum would be 
at the next higher order in t. The energies —/j, — t and 
— fj, + t are identical to the Taylor expansions of Eqs. J9|) 
and (fTDj) . respectively, to first order in t. The next order 
in the hopping amplitidc, second order perturbation the- 
ory would take into account hopping to an adjacent site 
with energy and back and would be of order t /fi, 
also what one would expect from expanding Eqs. §§§ and 
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FIG. 2: (color online) Density of states (DoS in arbitrary 
units, colored blue) versus energy E (in units of t) for a zig- 
zag domain wall a) and an armchair wall b) when /i = 0.5t. 
Also shown is the DoS for bulk bands (red). 



PU|) as well as (J5J) to second order in t. The order t 2 / [i 
contributions are momentum dependent and the energy 
levels become bands. Then t is made larger than /Lt, they 
spread out into the bands depicted in Fig. 2. 

In the neutral ground state, the half of the electron 
states with lowest energy will be filled. For the zig-zag 
wall there is a profound difference between two cases - 
when the midgap band Eq. (|10[) overlaps the negative 
energy bulk band (fi/t < 3/2) and when it doesn't (fi/t > 
3/2). When it doesn't, the neutral ground state has the 
negative energy bulk states Eq. (JS]) and the lower domain 
wall band Eq. © completely filled. There is a gap (which 
could be small) between the top of the filled lower bulk 
band and the empty upper domain wall band Eq. (I10p . 
When they do overlap (as depicted in Fig. 2), some of the 
upper domain wall band Eq. ((TO]) will fill before all of the 
negative energy bulk states are filled - the domain wall 
will borrow some electrons from the bulk. It will then have 
a finite charge density and a partially filled upper band 
which will behave like a one-dimensional metal, even in 
neutral graphene. 

The armchair domain wall is depicted in Fig. 1(c). It 
is oriented along the y-axis. It has the same gapped bulk 
branches Eq. j8]) as the zig-zag wall. We look for wave- 
functions which decay exponentially in distance \A X \, \B X \ 
from the wall. They are superpositions of two plane- 
waves propagating along the wall with wave-vectors k y 
and k y + 2tt/3. (This corresponds to mixing of the 
graphene valleys.) The spectrum has four bands depicted 

in Fig. 2. It is E = ±Jt 2 sin 2 ^k y + fi 2 K(k y ) where K ■ 
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k x and k x must be determined by solving two equa- 
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tions: K- 1 = £ (cosh ^fc^ 1 ' + cosh ^/ci 2) ) 2 and 



cosh^fsfci 2 ' 



COS ~2~ fey 5 



We 



can find explicit solutions in the large and small fi/t lim- 
its. When fi/t is large the spectrum is concentrated at 
four values 



fj,»t: E 



±(fX + t + 
±((l-t + 



(11) 



Two of these are inside the bulk spectrum and two are 
in the gap. They agree with what we would expect when 
(i is large, where there are two energy levels, —\x and 
+/i corresponding to electrons sitting on the black or 
white dots, respectively in Fig. 1(c). Then, the leading 
effect of turning on a small t is that the pairs of adjacent 
degenerate states that exist at the location of the domain 
wall are split by tunneling. Now, unlike for the zig-zag, 
there are degenerate pairs with both zeroth order energies 
+fi and — \x. The splitting produces four domain wall 
energies in Eq. (fTT|) . Further corrections are of order t 2 / n 
which, when taken into account spread the four levels into 
four bands, which then get wider as t gets larger. 
In the limit fx «t we also find four bands, 



fi «t : E 
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+ 4< 2 (1 + • 
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±Jt 2 sin 2 %k y 
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*2 (4_ COS 2 ifky) 2 



(12) 



As depicted in Fig. 2, the upper and lower band arc 
entirely within the upper and lower bulk bands. The 
middle two overlap the bulk energy gap and themselves 
have a gap which is much smaller than the bulk gap, 
AE = An 2 /it « 2fi. For a typical small /i ~ 25mev 
and t ~ 2.7et>, the gap in the mid-band states is tiny, less 



than lmev. This existence of a gap in the spectrum of 
states bound to the domain wall is compatible with the 
continuum analysis since the gap is vanishingly small in 
the continuum limit, scaling to zero with the lattice spac- 
ing a, so it is not visible to the continuum Dirac Hamil- 
tonian. Intuition for the gap in the armchair spectrum 
can also be gained by studying Fig. 1(c). Because of the 
alternating pattern of pairs of black and pairs of white 
dots as one follows the domain wall, the translation sym- 
metry along the domain wall is by two lattice spacings, 
rather than one. This reduced translation symmetry will 
gap the domain wall spectrum, analogous to gapping cre- 
ated by a Peierls instability. What is surprising here is 
that the gap is so small. 

In summary, we have shown that the simplest domain 
walls in gapped graphene can have interesting electronic 
properties. A partially filled domain wall band will be- 
have like a one-dimensional metal. The continuum anal- 
ysis suggests that similar behavior can be expected for 
other types of domain walls, such as those arising from 
reversing the pattern of a lattice distortion. Analysis 
of the details of the spectrum in those cases is left to 
future work. They could also occur in other materials 
which have a Dirac spectrum, such as the hypothetical 
flux phases of a square lattice where vortices have re- 
cently been discussed (26j and where domain walls, some 
of which would support zero modes should exist. 
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